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3  ABSTRACT 

A  frequently  occurring  problem  in  weapon  systems  analysis  is  the 
conpitation  of  expected  fractional  damage  of  an  area  target  engaged  by 
a  salvo  of  area  damaging  rounds.  Of  particular  interest  is  the  case 
Involving  both  round  to  round  and  occasion  to  occasion  errors.  When  the 
number  of  rounds  is  large,  the  available  solutions  may  encounter  acute 
computational  difficulty.  This  report  presents  a  computational 
procedure,  using  Jacobi  polynomials,  which  overcomes  this  difficulty. 
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a,b 


PX(K) 


PjOe.y) 

Pi4.(5>*n) 


LIST  OP  SYMBOLS 

Axes  of  elliptical  target  in  x  and  y  directions 
Expected  fraction  of  target  damaged 

Number  of  rounds  in  salvo 

Single  shot  probability  of  damaging  a  target  point  (x,y)  with 
a  round  impacting  at  (u^v^ 


Pp(u. v  )  Distribution  of  impact  points  of  the  i-th  round  about  its 
aim  point 


Distribution  of  (x,y) 

Distribution  of  occasion  to  occasion  errors 


"r 

<vV 

v(z) 

(x,y) 

(5/1) 


Target  radius  for  circular  target 

Impact  point  of  the  i-th  round 

Weighting  function  in  Jacobi  polynomials 
Target  point 

Cannon  aim  point  for  all  weapons 

Common  value  of  a  and  a  when  a  =  a 
x  y  x  y 

Parameter  in  the  damage  function,  p^(K) 


aVaT\ 

*r(z) 


Standard  deviations  of  the  round  to  round  errors 
Standard  deviation  of  the  occasion  to  occasion  errors 
The  r-th  Jacobi  polynomial 
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I.  INTRODUCTION 


A  frequently  occurring  problem  in  weapon  systems  analysis  is  the 

computation  of  the  expected  fraction  of  a  target  damaged  by  a  salvo  of 

area  damaging  rounds  all  aimed  at  the  same  aim  point.  The  many  variations 

of  this  problem  are  commonly  referred  to  as  coverage  problems.  Of 

particular  interest  is  the  case  involving  both  round  to  round  and 

occasion  to  occasion  errors.  The  case  involving  only  round  to  round 

errors  has  been  treated  by  Groves.^  The  more  general  case  involving 

both  types  of  errors,  in  the  form  considered  in  this  report,  is 

2 

contained  in  an  unpublished  work  by  Grubbs.  When  the  number  of  rounds 
is  large,  both  Groves’  and  Grubbs'  procedures  encounter  acute 
computational  difficulty.  An  analysis  of  this  difficulty  and  the 
presentation  of  an  alternate  computational  procedure  employing  Jacobi 
polynomials  comprise  the  substance  of  this  report. 

Let  p1(K)  =  p1(K| (ui,v1 ) ,  (x,y))  be  the  conditional  probability  of 
damaging  the  target  point  (x,y)  given  that  the  i-th  round  impacts  at 
point  (u^v^)  and  let  *>e  the  density  function  describing  the 

distribution  of  (u^v^)  about  the  aim  point.  If  all  rounds  in  a  valley 
are  identically  distributed  about  the  same  aim  point  ( 5,T|) ,  then  the 
probability  of  damage  over  all  impact  points  and  all  rounds  i,  i»l,2,...,N, 
is  equal  to  1  -  (l-z)^,  where 

OD  ® 

z  =  J  J  p^K)  P2(u,v)  dudv.  (1) 


If  one  assumes  (x,y)  is  distributed  over  the  target  area  as 
Tj(x,y)  and  the  intended  aim  point  (5,1))  is  itself  a  random  variable 
with  density  p^(5>"n)>  then  for  N  rounds  aimed  at  the  c  crane  n  aim  point, 
the  expected  fraction  of  the  target  damaged,  f  ,  is  giver,  by 


* 


References  are  lieted.  on  page  23. 
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(2) 


?N  *  J  J  J  J  [1  -  (1-z)N3  Pj(x>y)  P^(5,71)  dScLUdicdy 

X  «®  -0» 

where  T  is  the  target  area. 

The  difficulties  arising  in  the  solution  of  (2)  are  discussed  in 
Section  II.  A  computational  procedure  using  Jacobi  polynomials  to 
overcome  these  difficulties  is  developed  in  Section  III.  In  Section  IV 
an  illustrative  example  is  given  and  numerical  results  are  presented  in 
Section  V. 

II.  NATURE  OP  THE  PROBLEM 

The  procedure  used  by  Grubbs  and  Groves  to  solve  (2)  consists  of 
first  expanding  fN(z)  using  the  binomial  expansion 

t/z)  -  (l-z)»  .  £  <-l)J  (j)  ZJ.  (J) 

J=C 

Thus  (2)  can  be  written  as 


J-l 


where 

00  00 

Gj  =  J  J  J  J  P5(x,y)  Pj,(?,T1)  d^dTldxdy.  (5) 

For  some  target  distributions  and  for  some  forms  of  the  distributions 
P^(x,y)  and  p^ ( 5 , T( ) ,  G^  can  be  obtained  analytically,  see  e.g. 

Reference  [3];  however,  in  most  cases  approximations  are  required. 

In  some  problems  of  interest,  f  is  desired  for  large  N,  e.g., 

K  2  50.  Although  fy  lies  in  the  interval  0  s  s  1,  the  individual 
terms  of  the  series  in  (4)  can  become  extremely  large.  As  shown  in 
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U»\*tiun  V,  Table  1.,  the  partial  sumc  of  this  aeries  may  oscillate  in 
.'Ign,  Initially  with  increasing  magnitude  until  at  some  value  of  the 
.urination  Index  the  magnitude  begins  to  decrease.  For  large  N  and 
.-.  'n-  ,5.  the  quantity  j  Gj  may  attain  an  order  of  magnitude,  which 
ex.’cvds  the  word  length  of  a  digital  computer,  even  when  double  or 
triple  precision  -imputations  are  uced.  When  this  situation  is  present, 
the  roundoff  error  incurred  at  various  stages  of  the  computations  exceeds 
the  value  of  f.^,  hence  Equation  (k)  is  of  no  value  in  computing  the 
fractional  damage,  f.,. 

W 

The  source  of  the  computational  difficulty  is  the  binomial 
coefficients  arising  in  the  expansion  cf  (l-z)N.  Accordingly,  the 
approach  used  in  this  report  is  based  on  finding  an  alternate  expression 
for  (1-2 )N.  This  lc  done  in  terms  of  Jacobi  polynomials.  In  particular, 
seek  an  approximation  of  the  form 

M 

fH(a)  -  U-*)N  -  £  arN  »r(*),  0  *-  z  <  1,  M  *  N  (6) 

r«0 

where  $.(-)  io  the  r-th  Jacobi  polynomial  appropriate  to  the  weighting 
function 

w(z)  *>  z®(l-z)8  ,  cr  >  -  1,  0  >  -  1.  (T) 

Jlnce  (l-r.)5’  is  a  polynomial  of  degree  N,  Equation  (6)  is  an  exact 
representation  when  M  ®  r.'.  The  computation  of  f^(z)  using  Equation  (6) 
for  large  M  encounters  roundoff  problems  similar  to  those  present  w*-;>  n 
using  (3).  The  Jacobi  polynomials,  however,  provide  a  weighted  least 
squares  approximation  and  usually  lead  to  high  accuracy  ith  only  a  few 
terms.  This  is  particulary  true  for  the  Chebyshev  polynomials,  a 
special  case  of  the  Jacobi  polynomials.  These  polynomials  have  the 
advantage  that  each  successive  approximation  is  closer  to  f^Cz) ,  i.e., 

|Pr+l(z)  "  Vz)l  *  1Pr{z)  ’  fN(z)l  (8) 

for  r=0,l,.  ..,[1-1. 
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A  new  series  for  computing  employing  the  J&cobi  series  expansion 
of  (l-z)N  is  listed  in  Section  IV.  In  Section  V  the  behavior  of  this 
series  is  compared  to  the  series  derived  from  the  binomial  expansion. 
Equation  (4).  The  objective  of  the  studies  performed  was  to  determine 
whether  the  new  series  could  be  truncated  for  some  M  before  severe 
roundoff  error  arose,  and  still  retain  two  or  three  decimal  accuracy 

ln  V 


III.  EXPANSION  OF  (l-z)  IN  JACOBI  POLYNOMIALS 

The  Jacobi  polynomials  axe  orthogonal  in  the  interval  0  <  z  £  1,  the 
interval  of  irterest  for  z  as  defined  by  (l).  The  Jacobi  polynomial 
$r(z)  can  be  generated  by  the  expression 


«r(z)  =  cr  z^l-z)”9  U^(z), 


where 


uJ(*>  -  ^  [zrw(l-*rel, 
r  L  J 


and  c  is  an  arbitrary  or  normalizing  constant.  See,  e.g.,  Hildebrand, 
r 

The  coefficient  arN  in  (6)  is  obtained  from  the  integral 

arN  =  MyJ  J1  dz>  f11’ 

o 


where 


7r  =  J1  w(z)  i‘r(z) 


The  integral  in  (11)  can  be  integrated  by  parts.  Let  u  =  (l-z)”  and 
dv  =  Ur(z)  dz.  Then  du  =  -  N(l-z)^’^  and  v  =  Ur”^(z).  Since 


uj(0)  -  uj(l)  -  0  ,  j=l,2,...,r-l, 


I 


a 


rN 


(l-z)*’1  uj^U)  dz. 


After  integrating  by  parts  r  times  we  obtain 


(13) 


arN  =  (crAr)  N(N-l) . . .  (N-r+l)  f  U°(z)  (l-z)!‘'r  iz 

o 

=  [crN:/>r(N-r)’.]  J1  zr40f(l-z)‘J+S  dz 
o 

=  [crK’./rr(N-r)l]  B(r+flffl,  N+H+l),  (lU) 

1*. 

where  B(x,y)  is  the  beta  function-  As  stiown  by  Hildebrand  ,  page  271, 
the  quantity  y^  is  given  by 


7r  -  (-i)rr: 


r1 


rr 


>  J 


r+or,.,  Xr+S 
z  (1-z) 


dz 


*  (-l)rrl  Arr  cr  B(r-*or+l,  r-*fl+l),  (15) 

where  A  is  the  coefficient  of  zV  in  i  (z).  The  coefficients  A  .  in 
rr  r  rj 

i  (z)  can  be  obtained  by  use  of  the  hypergeoraetric  series.  See,  e.g., 

c 

Ccurant  and  Hilbert. J  $r(z)  is  given  by 


r 

Vz)  =  Y,  ArJ  zi>  <16) 

J=° 

where  A^^  is  obtained  by  replacing  s  by  (o+6+r+l),  t  by  -r  and  v  by 
( 1-tor)  in  the  j-th  coefficient,  Hj,  of  the  hypergeometric  series,  i.e., 

J 

H  =  n  (s+i-1)  (t+i-l)/(v+i-l)i,  J=l,2,...,r  (17) 

d  1=1 

and  H  =1.  Thus 

c 
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A  (-lpr-^+e  +  p(r*Or+6+2)...(2r-*o+6'  / 

Vr  (l-io)(2+o)...(r-+o)  < 

(-l)rr(2r+o+6+l)r(l+o'  , 

r(r+o+0+l)  ' 

whor<-  P(x)  io  th*  iJamma  function.  Inserting  the  results  of  (15)  and 
fl?)  Into  (iU)  ani  employing  the  relationship 


we  obtain 


B(x,y>  r(x)r(y)/r(x-*y) 


( N  (2r+o+0+l)r(r-«a+0+l)r(r-«o+l5r(N+e+l] 


rN  \rt  r(i-»a)r(N+r-K>+B+2)r(r+6+l) 

Equation  (20 )  Is  made  more  amenable  to  computation  by  noting  that 

(8+l)(B+2)...(B+N) _ 


oN  (0+6+2) (0+6+3) • • • (o+B+N+lJ 


•hi  *  »']•  <22> 

Tv  avoid  an  indetermi  nancy  for  r*l,  o  +  8**l,  a^  should  be  written 

•»  -  ■^CTis}(arfi)]- 


IV.  ILLUSTRATIVE  EXAMPLE 


b  t  T  -  |(x,y)|xc/a‘  +  y2 /b2  r-  l}  , 

P1(K|(u1,vi),  (x,y) )  =  exp  •[(-l/2o2)[(x-u1)2  +  (y-v*)2]}  >  (23) 

P2(ui»vi)  =  (l/2noxo  )  exp  •[(-l/2)[(ui-5)2/<*2  +  (v1“-n)2/<7y]}  >  (2l+) 


P5(x,y)  =  1/nab 


14 


Pi/s.l)  *  (l/^c^o^)  exp  {(-l/2)[(5/o^)2  ♦  ("n/c-ti)2]}  *  (26) 

For  the  above  example  the  common  aiming  point  (?,Tl)  is  distributed 
about  the  target  center,  since  the  target  center  is  (0,0)  and  the  mean 
of  is  (0,0).  By  inserting  the  above  representations  into 

Equation  (l)  it  can  be  shown  that  z  reduces  to 

z  =  q' exp  {(-l/2)[(5-x)2/(c2+c2)  +  (‘n-y)2/(c2+c2)]}  ,  (2?) 


where 


’  '  vC(ikm«)(<’iK)]  • 

From  (2),  (5),  (6)  and  (l£) 


W-I 


L.  rN  L 

r=0  J=0 


A  .  G, 

rj  J 


where  G  is  given  by  Equation  (5).  By  inserting  the  distributions 
J 

specified  by  Equations  (25)  through  (28)  into  (5),  by  completing  the 
squares  in  z^p^(f,T'l)  and  integrating  over  5  and  "H «  G^  can  be  shown  to 
reduce  to 

Gj  -  C(qj/J)  IjJ  j=l,2,...N,  ( 


•-he  re 


=  -*■> 

=  44+0x)(cK+cy}]  '  /ab’ 

(31) 

=  [ ( J  a^+o2+o2 )  /J  J  7  , 

(32) 

=  [(j°^VCTy)/j]  7  ’ 

(33) 
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Ij  -  (1/2*!^  vJJ  exp  [(-l/2)(x2/D^y2/Ej)]  dxdy,  (3M 

T 

T  being  the  ellipse  x2/a2  +  y2/b2  =1. 

Since  A  *  1,  Equation  (29)  can  be  written 

N  Nr 

?»  '  1  -  I  *r»  '  C  I  *r«  I  V  <«> 

r=0  r=l  J=1 

2 

The  aeries  developed  by  Grubbs  ,  using  the  binomial  expansion,  takes 
the  form 

N 

VC  I  (*DJ+1(j)  CqJ/j>  Xj •  (3 6) 

J=1 


When  *  0,  i.e.,  no  target  location  error,  this  series  is 

identical  to  that  described  by  Groves^ . 

Analytical  approximations  for  I,  are  contained  in  the  paper  by 
2  " 

Grubbs  .  I,  can  also  be  obtained  by  numerical  methods,  e.g.,  by  the 

J  g 

method  described  by  Breaux  or  alternatively  as  follows:  Let 


erf 


(57) 


On  many  computers  this  function  is  as  standard  as  the  trigonometric  or 
exponential  functions  and  hence  can  be  used  to  eliminate  one  integral 
in  (3M-  Equation  (3*0  can  then  be  written 

Ij  -  f  exp  [-(l/2)y2/E2]  erf  [a(l-y2/b2)l/2/(V?  Dj)]  dy.  (38) 


Ij  is  now  in  the  form  of  a  single  integral  and  can  easily  be  obtained 
on  a  digital  computer  by  use  of  standard  subroutines. 
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Whf  ii  ^  -  f  ~  and  tor  circular  suir,<  ts,  I.*  ., 


t  -  1<y»  w*  obtain  the-  familiar  result 


1  .  =  [l  -  Xp  (-R^/rV)]  , 


.‘l'  r-  o  ir.  the  common  value  of  D,  and  E,.  When  o  -  0,  I,  taker. 

J  J  J  5  j 

th*.*  form 

Ij  J  U  -  '*xp  (-'J)]  ( 


where 


X  -  !^/£(o^c2) 


To  provide  an  independent  method  for  •‘.eeking  the  accuracy  and 

convergence  of  the  ooriev.  developed  in  thin  report,  the  numerical  .studio;. 

wer<  performed  for  euros  whor*  I.  .‘ou  Ld  b»-  v-  pr* .  <nt<.*u  by  Equation  (4»  ). 

1  **  1 

For  thio  ease  Grubbs  sorle~  is  identical  t  Grove a  and  can  be  written 

N 

7u  “  c  I  ("i)J+1  (j)  (Ul) 

J=i 

7 

For  thio  cuce  Breaux  has  found  an  alternate  series 


fN  **  c  Y.  [(^’V  -  U~q)J]/J-  (**« 

J=i 

Thin  series  is  eaoily  summable  for  all  valuer  of  the  parameters,  av,  a, 

R^,,  and  N  and  provides  an  exact  solution  for  comparison.  Note  that  the 

above  simplifications  ore  made  only  to  provide  an  exact  solution  for 

numerical  comparison.  The  general  structure  of  the  computational 

procedure  does  not  necessitate  that  I,  be  representable  by  (40)  since 

J 

as  pointed  out  previously  the  integral  form  of  Ij  in  Equation  (38)  can 
be  obtained  in  fractions  of  a  second  by  use  of  standard  techniques  for 
numerical  quadrature. 
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V.  NUMERICAL  RESULTS 

The  numerical  studies  were  designed  to  examine  the  accuracy  and 
rate  of  convergence  of  the  following  series: 

4X)  =  C  £  (-1)J+1  (jj  (qJ/j)  Ij  (U3) 

j=l 

M  Mr 

42’  ■  1  -  1  arN  -  C  1  “r»  I  *J  <“> 

r=0  r=l  J=1 

M  r 

'  "  43>  -  -  C  1  arH  I  rj  <k5) 

r=l  J=1 

Equation  (U3)  is  the  series  arising  from  (Ul),  (hk)  is  that  arising  from 
(35) •  By  setting  z  -  0  and  M  =  N  in  (6)  it  can  be  shown  that 


I  arN  =  lm 


(3) 


Applying  this  result  to  Equation  (35)  leads  to  the  third  series 

Numerical  experiments,  see  Table  I,  indicate  that  when  f  is  not 


s(2) 


(3) 


equal  to  0  or  1.0,  decreases  monotonically  and  SjJj  '  increases 

monotonically  as  M  increases,  both  approaching  f^.  A  weighted  average 


?(2) 


s(5) 


of  the  two  sequences  and  '  would  thus  seem  to  offer  a  more 

accurate  approximation  to  f^.  Such  a  new  sequence  can  be  constructed 
as  follows :  Let 


*4°  - I41’  -  41-i|  ■ 


Then  define 


410  ■  [42)/i42)  +  43>/a43)]/[i/'i42)  +  i/a45)]- 
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fin 


weight. 
!  \ 

!  n  • 


* 

-IK  ? 


uttu  *h- J  to  -  u  *h  sequence  are  Inversely  proportions:  to 

th  o  .in*  irt  with  smaller  would  seemingly  b*.-  tl*-  su-t 

/ 1  ^  M 

AlterntiUv*  iy,  .1,  oun  be  view  d  a.'  the  qu«  rv-"  farmed  by 

‘4  „(?) 

*tl  .r  th*-  linear  extrapolation  of  i  he  tv<>  s-.quen  •••.:  3,^ 


K«  ;  ulti  of  numerical  exp*  riments,  not  listed  her*.  ,  indicate  that 
hf.  t.  convergence  Is  attained  for  weight  functions  having  a  -  B  =•  -  l/2, 


w( u )  »  z"  l^2  (l-z)‘  1^* 


For  this  weight  function  the  Jacobi  polynomials  reduce  to  the 
"shifted"  Chebyshev  polynomials.  These  polynomials  generally  have  th*. 
br.it  convergence  properties,  as  verified  in  this  case  by  experiment. 


Discussion  of  Results 

The  critical  parameters  effecting  the  accuracy  and  convergence  of 


the  series  solutions  are  a  and  N.  When  *=  0,  and  a x 

Equation  (28)  can  be  written 


=  o. 


o“  (l-q)/q. 


(U7) 


The  cases  studied  were  for  *  i.0,  q  =  .01,  .1,  .5,  -5>  *7,  .9/ 

::  -  50,  luO,  15b,  200,  500,  Uoo,  500,  1000,  and  TR  -  50  with  c  constrained 
by  Equation  (1*7). 

A  comparison  of  the  convergence  of  the  four  series,  Equations  (U$) 

through  (U6)  is  illustrated  for  a  typical  case  in  Table  I.  The  series 

arising  from  the  binomial  expansion  is  seen  to  oscillate  with  extremely 

large  magnitude  and  does  not  provide  a  useful  result  until  most  of  the 

50  terms  are  added.  S,[2^  and  on  the  other  hand,  approach  the 

true  solution  monotonlcaliy  as  does  ,  the  weighted  average.  It 

should  be  noted  that  the  computations  were  performed  on  BP.LESC 

(Ballistic  Research  Laboratories  Electronic  Scientific  Computer)  which 

56  IT 

has  a  useable  word  length  equivalent  of  2?  or  approximately  10  .  The 

normal  word  length  on  most  commercial  machines  is  approximately  10 
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For  this  reason  It  seems  unlikely  that  could  be  summed  at  all, 

for  this  case,  on  a  connercial  machine,  except  by  use  of  double 
precision. 

Table  II  is  a  listing  of  the  parametric  study  to  compare  the 
accuracy  of  the  four  series.  For  each  pair  of  values  ( II , q ) ,  the  five 
entries  in  succession  are  the  exact  solution,  S°  followed  by 

SA  ,  i=l, 2,5,U.  The  asteriks  (*)  denote  that  no  useful  result  could  be 

™  (]  )  (i) 

obtained  from  .  The  sequences  ,  1=1, 2, 5,^  were  terminated  either 

when  M  =  N,  or  when  any  intermediate  number  exceeded  10^.  5y  inspection 

of  Table  I  it  is  seen  that  is  the  most  accurate  scries  with  a 

maximum  error  of  .0002.  Note  that  31  of  the  48  entries  could  not  be 

computed  using  the  series  derived  from  the  binomial  expansion.  The 

average  computation  time  for  the  cases  studied  was  approximately  1.1+ 

seconds  per  case  on  BFLESC. 

VI.  CONCLUSIONS 

A  computational  procedure  for  determining  expected  fractional 
damage  for  an  area  target  engaged  by  a  salvo  of  area  kill  weapons  has 
been  presented.  The  procedure  employs  Jacobi  polynomials  and  in  most 
cases  the  successive  approximations  converge  rapidly  to  the  true 
solution.  Two  new  series  solutions  have  been  presented,  one  increasing 
monotonically  and  the  other  decreasing  monotonically  with  the  summation 
index,  both  approaching  the  true  solution.  A  method  for  averaging  the 
two  solutions  has  also  been  presented  which  accelerates  the  convergence, 
thereby  making  the  method  useful  even  in  extreme  cases  where  numerical 
difficulties  force  the  termination  of  the  series  before  convergence 
of  either  has  been  reached. 
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* 

.40000000-  1 

.92041076  * 
•76460313 

* 

.24969175-  3 

* 

.14726577-  2 

-.401000k. 

.61806674 

.76752166-  3 

.29402156-  2 

•  38326<ao  1 

.40542272 

.13043077-  2 

.32562751-  2 

-.2974514..  2 

.37004504 

.18414197-  2 

.354751-7-  2 

.lJ2673f'l  5 

.27358342 

.23304005-  2 

. 56965226-  2 

-.10584329  4 

.19609150 

.27434043-  2 

.57690745-  2 

.50081576  4 

.1362  -484 

.30847220-  2 

.56367494-  2 

-.20670582  5 

.91970193-  1 

.33441125-  2 

.36607293- 

.75196713  5 

.60424990-  1 

.35370562-2 

.3=626 a 9-  2 

-.24317858  6 

.36873154-  1 

.36729533-  2 

.369552C6-  2 

.7040536 )  6 

.24742770-  1 

.37640344-  2 

1  .5£?85923-  2 

-.18557128  7 

.I5656IC7-  1 

.36231318-  2 

.3>02625S-  2 

.43318622  7 

.10496434-  1 

.36593152-  2 

.39038179-  2 

-.92902972  7 

.73983605-  2 

•368C7iee-  2 

.39C45565-  2 

. 18171976  8 

.56827172-  2 

.56927868-  2 

-  3S<552695-  2 

-.32515228  6 

.47729915-  2 

.38902717-  2 

.39054557-  2 

-53354859  8 

.43114050-  2 

.39026279-  2 

.39055786-  2 

-.80450360  8 

.40874159-  2 

.39042727-  2 

•39C 5607S-  2 

.11169954  9 

•39835488-  2 

.39050474-  2 

.39056255-  2 

- . 14295696  9 

-  393755'-fi-  2 

-39053948-  2 

.59036358-  2 

.16886754  9 

.39161221-  2 

.35055429-  2 

O'90 5c;cl-  2 

-.18425666  9 

. 39IO296O-  2  . 

,  -39056036-  2 

.39036397-  2 

.10582284  9 

1 

-.17326003  9 

a  *  9  *  -  .5 

•  1U?-396U7 
-.11909445 
.87743005 
-.59707817 
•57494008 
-.2170  1  903 
.11561148 
-.565632tl 
.25387117 
-.10414660 
•58943608 
-.132246  >6 
.4C605&00 
-.11214117 
.27680621 
-.60601240 
. 11657346 
-.19464392 
.27843075 
-.52800155 
.36381647- 
.14242527- 
.40448705- 
.39005318- 
.59057258- 
6342- 
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0,  =  ° 


Exact  Solution  =  .39056596-02 


This  column  of  numbers  indicates  the  power 
of  ten  by  which  to  multiply  the  table  entry 

to  obtain  S^ . 


TABLE  II 

PARAMETRIC  STUDY  OF  SERIES  ACCURACY 
Of  a  8  *  -  l/2,  0  £  -  Gji  a  0 ,  Rrj,*  50,  0^  =  1 


.01 

.1 

.3 

•  5 

.7 

C 

*  y' 

s° 

.'0355B 

.01756 

.00879 

.00609 

.'00475 

.00391 

s1 

.05558 

.01758 

.GC&79 

.OOfCO 

.00473 

.00391 

50 

s3 

.03558 

.01756 

.00879 

•OC6C9 

.00473 

.00391 

s3 

.03558 

.01758 

.00879 

.OObOQ 

.00473 

.00351 

s4 

.03558 

.01758 

.00879 

.OC609 

.00473 

.003/1 

s° 

.06 38 3 

.02508 

.01062 

.0071- 

.00 5 52 

.00452 

s1 

.06385 

.02508 

.01062 

* 

* 

♦ 

100 

s3 

.06385 

.02508 

.01062 

.00724 

.00602 

.005 >i 

s3 

.06363 

.02308 

.01062 

.00719 

CO1'  0,0 

.00451 

s4 

.06383 

.02308 

.01062 

.00719 

.0055 2 

.00452 

s° 

.08673 

.02633 

.01170 

.00784 

.0055-8 

.00486 

s1 

.08673 

.0265- 

♦ 

* 

* 

* 

150 

s3 

.CS873 

.026 *1 

.01174 

.00879 

.OK,53 

.0l4c6 

s3 

.08673 

.02^31 

.01170 

.00782 

.00594 

.00480 

s4 

.08673 

.02631 

.01170 

.00784 

.00598 

.00487 

s° 

.10566 

.02860 

.01246 

.00829 

.cc6 51 

.005 15 

s1 

. 10566 

.02860 

* 

* 

*  ■ 

* 

200 

s3 

. 10566 

.02860 

.01283 

. .01255 

.02035 

.02927 

s3 

.10566 

.02860 

.01246 

.00824 

.00617 

.0C4O4 

s4 

.10566 

.02860 

.01246 

.00829 

.00631 

.00513 

s° 

.13522 

.03184 

.01354 

.00894 

.00677 

.CC540 

s1 

.13522 

* 

* 

♦ 

* 

* 

300 

s2 

.13522 

.03184 

.01717 

.02855 

.05168 

.07106 

s3 

.13522 

.03184 

.01346 

.00866 

.00651 

.00496 

s4 

.13522 

.03184 

.01554 

.00893 

.00675 

O 

O 

vn 

£- 

-4 

s° 

.157L7  . 

.03414 

■ .01431 

.00540 

.00711 

.00575 

s1 

.15747 

* 

* 

* 

W 

* 

LOO 

s2 

.15747 

.03416 

.02605 

.05259 

.08938 

.11640 

s3 

.15747 

.C3414 

.01403 

.00678 

.00625 

.00484 

s* 

.15747 

.03414 

.01450 

.go  938 

. CC7 06 

.00570 

s° 

.17510 

.03592 

.01490 

.  00  97 6 

.00736 

.00594 

s1 

.17510 

* 

* 

* 

* 

* 

500 

0 

.17510 

.03604 

.03903 

.08015 

.12731 

.15958 

s3 

.17510 

•03591 

.01433 

.00674 

.cc  6 10 

.C0468 

s4 

.17510 

.05532 

.01486 

.00971 

.00729 

.00587 

s° 

.23042 

.04146 

.01675 

.01087 

.00815 

.CO656 

s1 

.23042 

* 

* 

* 

* 

* 

1000 

s3 

.25119 

.04780 

.12643 

.21036 

.27647 

.31823 

s3 

.22958 

.04099 

.Cl4o8 

.00791 

.00525  v 

.00893 

s4 

.23039 

.04145 

.01661 

.01067 

.00795  ; 

.00635 
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